THE EXPERIMENTAL DESIGN AND OPERATION 
OF A ROTATING WICKLESS HEAT PIPE 



by 



Thomas James Daley 



LIBRARY 

NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIF. 93940 



United States 
Naval Postgraduate School 




THESIS 



THE EXPERIMENTAL DESIGN AND OPERATION 
OF A ROTATING WICKLESS HEAT PIPE 



by 



Thomas James Daley 



June 1970 



Tfvu document ha& been approved (JoA pabZia fie- 
tejxbe and &aZe; Zt6 dl6tAibvution untmited. 



i 



THE EXPERIMENTAL DESIGN AND OPERATION 



OF A ROTATING WICKLESS HEAT PIPE 



by 



Thomas James^Daley 

Lieutenant, Junior Grade, United States Navy 
B.S., United States Naval Academy, 19^9 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OP SCIENCE IN MECHANICAL ENGINEERING 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1970 



OSTGRADUATE SCHOOD 
:]Y, CALIF. 93940 



ABSTRACT 

An experimental rotating wickless heat pipe apparatus was 
designed and machined. The apparatus includes a rotating heat pipe 
assembly, test stand, spray cooling assembly, safety shielding, 
and instrumentation. 

A revised condensing limit for the operation of the rotating 
heat pipe was obtained by modifying Ballback's Nusselt film con- 
densation theory to Include the effects of a thermal resistance in 
the condenser wall and in the condenser outside surface cooling 
mechanism. Approximate results, obtained for half-cone angles of 
1, 2, and 3 degrees, show that less heat can be removed than originally 
predicted by Ballback[7 ] , and that the outside heat transfer co- 
efficient can significantly alter the condensing limit. 

An improved Nusselt theory was developed which applies for 
all half-cone angles, and which includes the effects of the thermal 
resistances in the condenser wall and in the condenser outside surface 
cooling mechanism. This formulation led to a second-order non-linear 
differential equation for the film thickness which was numerically 
integrated using a free-overfall boundary condition at the condenser 
exit. Results obtained for a half- cone angle of 0 degrees are sub- 
stantially less than the results obtained from the approximate solution 
for half-cone angles of 1, 2, and 3 degrees. 
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I , INTRODUCTION 



A, CONVENTIONAL HEAT PIPE 

A heat pipe is a self-contained device capable of transporting 
large quantities of heat at nearly isothermal conditions. Con- 
ventionally it consists of three main parts: a container, a working 

fluid, and a wicking structure capable of promoting capillary action. 

When heat is added to one end of the container, some of the 
working fluid evaporates. This end is the evaporator. The pressure 
differential in the vapor chamber due to the increase in vapor at 
the evaporator causes the vapor to flow to the opposite end of the 
container, transporting heat as latent heat of vaporization. Here 
the vapor condenses on the cooler surface and releases the heat to 
be removed from the structure. This end of the container becomes 
the condenser. The condensate is pumped back to the evaporator 
section by capillary action, thus completing the cycle. The heat pipe 
is unique in that it operates with no moving parts, maintains nearly 
Isothermal conditions, and removes the dependency on gravity through 
the use of capillary pumping. 

Limits on the heat transfer capabilities of conventional heat 
pipes are imposed by a number of fluid dynamic mechanisms , including a 
wick resistance limit, a sonic vapor velocity limit, a vapor velocity 
entrainment limit, and also ? wick boiling limit. These limits have 
been considered in a number of studies [1, 2, 3, 4 ]* with the con- 
clusion that in general a conventional heat pipe using an ordinary 
working fluid is limited by the rate that the condensate can be returned 
to the evaporator due to wick resistance, 

* Numbers in brackets indicate references listed in Bibliography, 
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B. THE ROTATING, WICKLESS HEAT PIPE 

The restriction imposed by the inability of the wicking 
structure to provide sufficient pumping force to return condensate 
to the evaporator can be alleviated by removing the wick and 
utilizing centrifugal acceleration as a pumping source • If the 
pipe is allowed to rotate about its longitudinal axis and the wall 
is given a slight taper, a rotational force field can easily be 
produced that not only overcomes the viscous losses in the system, 
but also the losses due to a hindering gravitational field. The 
operation of a rotating heat pipe is shown schematically in FIGURE 1. 




FIGURE 1. Schematic Drawing of Typical Rotating Heat Pipe 

The applications of a rotating heat pipe appear to be wide and 
varied. Wherever there is rotation and a need to transfer heat, the 
use of such a system is feasible. The cooling of motor rotors, turbine 
rotor blades, aircraft braking systems, and air-conditioning units are 
just a few examples of the uses to which the rotating heat pipe concept 
can be applied. 
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The capillarity of the wi eking structure need no longer be 
considered as a limit on heat pipe operation since it has been 
replaced as a pumping force by an artificial gravity head, Ballback 
[7] has studied the limitations imposed on the rotating system by 
various fluid dynamic mechanisms in a preliminary theoretical analysis 
of the operation of a wickless heat pipe. By using existing theories 
and experimental correlations, he was able to estimate the limits on 
rotating heat pipe operation imposed by the onset of critical nucleate 
boiling, condensate entrainment, and choked vapor flow. An estimate 
of a condensing limit was obtained by performing a Nusselt-type analysis 
for condensation on a truncated cone, since existing theory and 
correlations were not applicable. This analysis, however, was restricted 
by certain simplifying assumptions that were made to obtain a closed- 
form analytical result, Ballback ’s results indicated that the operation 
of a rotating heat pipe is limited by the amount of heat that can be 
removed from the condenser section, 

C. THESIS OBJECTIVES 

The objectives of this study were to: a.) design and construct 

an experimental rotating heat pipe assembly; b,) revise and supplement 
the initial theoretical study of rotating heat pipe operation performed 
by Ballback, 
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II. EXPERIMENTAL PROGRAM 



A. BASIC DESIGN CONCEPTS 

For Ballback to generate numerical results from his theoretical 
model, a fixed geometry for the heat pipe had to be decided upon. 

He chose to model the evaporator section along the lines of the rotating 
boiler used at the Lewis Research Center in Cleveland, Ohio, for the 
study of boiling heat- transfer coefficients operating under high gravity 
levels • With this geometry fixed he allowed the size of his 
condenser section (length and base radius) to be governed by the amount 
of taper in the section. The taper was defined by the half-cone angle, 
which was a parameter in the theoretical study. A strong influence in 
the overall design was the desire to be able to visually observe the 
mechanisms that were taking place within the heat pipe during operation. 
The onset of critical nucleate boiling and the entrainment of the 
returning condensate were two principal phenomena which were to be 
observed . 

B. DESCRIPTION OP EQUIPMENT 

The components of the wickless heat pipe are divided into four 
sub-groups to simplify the explanation of their design and construction. 
These are: 1.) evaporator 2.) condenser, 3.) auxiliary equipment, and 

4.) instrumentation. FIGURE 2 is a schematic representation of all the 
components of the test apparatus. A cross-sectional drawing of the 
assembled heat pipe is pictured in FIGURE 3 and FIGURE 4 is a photograph 
of the machined pieces of the rotating heat pipe prior to assembly. 
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FIGURE 2 SCHEMATIC DIAGRAM OF APPARATUS 
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FIGURE 3 CROSS SECTION OF ROTATING HEAT PIPE 
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1 . ) The Evaporator 

The evaporator was modeled to approximate the rotating 
boiler tested at Lewis Research Center [6] . It was machined of 
stainless-steel because of this metal's corrosion- resistant properties. 
It has an Inside diameter of 3-125-lnches, an exit diameter of 2.625- 
Inches, and an overall length of 5-905-inches (see FIGURE 3). This 
length Includes 2/3-lnch for ‘the Insertion of two pyrex-glass view- 
ports, 4-lnches of heating surface, and an Intended -g-lnch adiabatic 
section. The Inside diameter of the evaporator Is stepped |--lnch from 
the exit diameter to allow for the formation of a liquid annulus. The 
Inner viewport rests on a teflon-coated metallic 0-rlng seal and Is 
separated from the outer viewport by a 3/32-lnch thick compressed-fiber 
gasket. The viewports are held In place by a stainless-steel end cap 
which Is secured In place by twelve v-inch stud-bolts. These bolts 
also support two bronze collector rings. The rings are Insulated from 
their supports by l/l6-lnch teflon tubing and washers. Electrical power 
is passed by a brush assembly consisting of eight spring-loaded graphite 
brushes through the collector rings to an 11-gage Chromel-A heater wire 
in a 3/l6-inch inconel sheath. Power to the heater Is furnished by a 
DC motor-generator capable of delivering I50 amperes at 220 volts. The 
heater was helically wound and silver-soldered In place to insure good 
thermal contact. After the heater was coated with a thin layer of 
Sauereisen cement, asbestos Insulation was packed around the outside 
to reduce the radial heat losses. In addition two radial l/l6-lnch 
grooves were machined on the outside surface to within l/l6-lnch of 
the Inside boiler wall at either end of the heating element to lessen 
axial heat losses. 
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With this construction and by assuming a maximum operating 

temperature for the heater element, it is possible to estimate the 

maximum heat flux into the boiler. Likewise, it is possible to 

estimate the heat losses by approximating the outside convection 

coefficient from an analysis of air flow over a stationary cylinder 

[11 ] . These calculations are shown in Appendix A. The results 

indicate that a maximum radial heat flux of approximately 470,000 
2 

Btu/hr-ft can be obtained with an associated radial heat loss of 

2 

approximately 1300 Btu/hr-ft (with 3/4-inch asbestos insulation) . 

As noted in Appendix A, the heat flux is dependent upon the 
thermal resistances that lie in the radial heat transfer path. The 
dominant resistance is that of the stainless-steel evaporator wall 
because of its low thermal conductivity. It is therefore advantageous 
to keep the thickness of the wall as small as possible. However, as 
mentioned in the section on Instrumentation, thermocouples are to be 
imbedded in the boiler wall to monitor the temperature to obtain a 
temperature gradient and associated heat flux. Prom an uncertainty 
analysis (Appendix B) , it is found that the farther apart the thermocouples 
are radially mounted to one another, the smaller the uncertainty in 
determining an accurate heat flux. So there is a trade-off on the 
selection of the wall thickness. In this case, a thickness of ^-inch 
was used. 

A photograph of the machined evaporator section is shown In 
FIGURE 5- 
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FIGURE 5 Machined Evaporator Section 
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2 . ) The Condenser 



The condenser section was machined of stainless-steel in 
the shape of a truncated cone with a 3 degree half-cone angle. 

The base of the condenser mates to the evaporator exit and has an 
inside diameter of 2.5-inches (see FIGURE 3). A step is created 
at the evaporator— condenser intersection where the condensate flows 
from the condenser to the recessed inside diameter of the evaporator. 

The condenser tapers from an inside diameter of 2.5-inches to one of 
1.46-inches at the truncated end over a length of 10-inches. From here 
the condenser extends an additional 1 1/2-inches to house a cylindrical 
end plug. Of the overall length, 9-inches (along the axis) is con- 
sidered to be the cooling surface. The first inch of the condenser 
surface near the evaporator is assumed to be adiabatic. The wall 
thickness of the condenser is 1/ 16-inch to keep the thermal resistance 
as small as possible. 

The condenser section has two flanges. The large front flange 
mates to the bearing-scat flange on the evaporator. They are held 
together by the previously-mentioned 1/4-inch stud-bolts. A flange at 
the small end mates to a similar flange on the end plug, through which 
there are eight 1/8- inch stud-bolts for assembly purposes. The end 
plug has an extention on its open end that has four 1-inch square 
surfaces upon which are mounted four phenolic junction boards. 

A photograph of the condenser and end plug (with shaft in place) 
is shown in FIGURE 6. 

3. ) Atixiliary Equipment 

The auxiliary equipment are sub-divided into four groups 
for ease of description. These are: a.) drive assembly, b.) test 

stand, c.) spray cooling assembly, and d.) safety shields. 
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FIGURE 6 Machined Condenser and End Plug (and Slip Ring Unit) 



a. Drive Assembly 

The drive assembly consists of the end-plug, shaft, 
pulley, drive motor, and bearing. The end plug has a 3/4- inch 
threaded shaft way machined into it (See FIGURE 3)- Between the 

shaftway and end of the end plug is a .27 inch hole to allow for the 

escape of non-condensable gases during filling (to be described later) 
and for the insertion of a pressure transducer into the vapor space 
during operation. The shaft is a stainless-steel cylinder which is 
hollowed to allow the passage of thermocouple leads from the junction 

boards to a slip-ring unit. The shaft has 3/4- inch of its outside 

diameter threaded for screwing into the shaftway, and is threaded for 
3/8-inch along its inside diameter on both ends. The condenser end 
holds a ^-inch arm that extends through the hole in the end plug into 
the vapor space, and into which is screwed a small, semi-conductor pressure 
transducer, using a teflon-covered thread for sealing. The opposite 
end of the shaft is threaded to accomodate a coupling to the slip-ring 

unit that rotates with the pipe. 

One-half of a square keyway is machined on the shaft and 

one-half on the shaftway such that they coincide when the shaft is 
screwed securely into place in the end plug. A stainless-steel square 
key is slid into the keyway. The key transmits the rotational torque 
from the shaft to the heat pipe. The torque is applied to the shaft 
via a 2 . 65 -lnch pulley and v-belt by a 2 -horsepower, 3 -phase variable 
speed motor. The motor has a maximum speed of 4500 rpm and a minimum 
speed of 450 rpm. The speed is controlled by a push-button electric 
remote control unit. Ihe motor is equipped with a magnetic disc brake 
for rapid shut-down. 
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A double rov;, angular contact bearing is mounted on a 
seat machined on the shaft. This bearing and the large .^6-lnch 
inside diameter), floating, single bearing (with bakelite compound 
cage to withstand high rpm) mounted on the bearing seat on the evap- 
orator flange allow the pipe to easily rotate in any orientation, 
b. Test Stand 

The test stand was designed to sturdily support the 
heat pipe assembly mounted on a 15x44-inch steel bed-plate and the 
motor on its bed-plate. The system operates in orientations ranging 
from 0 to 90 degrees from the horizontal . Both bed-plates are welded 
to a 4-inch iron pipe. The pipe is placed in three- 2-inch thick 
clamps which are supported 2-feet off the ground by v-inch steel plate 
As the pipe is rotated to any orientation, the motor rotates with it 
and the assembly is then clamped securely in place. Vibrations are 
kept to a minimum by vibration-mounting the motor. FIGURE 7 is a 
drawing of a side and front view of the test stand configuration. 



Heat Pipe Motor 




Side Front 

FIGURE 7 Schematic Drawing of Test Stand Configuration 
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c. Spray Cooling Assembly 

The cooling of the condenser section is accomplished 
by spray cooling. The assembly consists of a cylinder rolled from 
1/8-inch thick stainless-steel with an outside diameter of 13-inches 
that encloses the condenser section. The cylinder is cut in half 
longitudinally. The bottom half is secured to the bed-plate on which 
the heat pipe is mounted and the other half is bolted to its mate 
once the heat pipe is in place and ready for operation. Each half cylinder 
holds four spray nozzles at the same axial position, ihe spray nozzles 
are easily interchanged so that mists of droplets ranging in size 
from 300 to 600 microns can be obtained. To keep the spray from 
interfering with the operation of the heater or either of the bearings, 
1/8-inch thick shields are welded to the ends of the half cylinders. 

When the cylinders are mated these shields come within l/l6-inch of 
the bottom and sides of ridges machined on both the front and rear 
condenser flanges (see FIGURE 3). 

The coolant flows from copper feed lines to a steel mixing tube 
mounted on top of the spray cylinder. Here the average temperature 
of the coolant is measured by a quartz thermometer. Prom the mixing 
tube the coolant flows through plastic tubing to each of the nozzles. 

Once sprayed it passes through a drain in the cylinder and bed-plate 
to a second mixing tube which is also monitored by a quartz thermometer. 

Prom here the coolant is sent to an external drain, 
d.) Safety Shields 

The entire heat pipe is surrounded by l/8-lnch 
stainless-steel shielding to insure safe operation. The shielding consists 
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of three components: the evaporator housing, condenser housing, 

and shaft housing. 

The evaporator housing Is an Inverted-U shaped section 
that encloses the evaporator. The front end Is closed and contains 
a plexl-glass window to observe Internal pipe operation through 
the pipe viewports. The back end Is open and mates with the spray 
cooling cylinder end shields. The siding has a curved plexiglass 
window for the visual observance of exterior pipe operation and a 
door located near the bed-plate for the adjustment of the brush 
assembly. The Inverted-u section slides onto set screws on the bed- 
plate for easy mounting. 

The condenser spray cooling cylinder serves not only as 
the spray cooling assembly but also as the condenser housing safety 
shield . 

The shaft housing Is also shaped like an Inverted U and 
encloses the shaft, pulley, and rear bearing assemblies. Both ends 
are left open with one sliding against the spray cooling cylinder 
end shields. A slot Is left In the siding for the operation of the 
v-belt. This section also slides onto set screws on the bed-plate. 

4. Instrumentation 

The Instrumentation Is grouped Into the following sub- 
divisions for ease of explanation; a.) temperature measurements, 
b.) pressure measurements, c.) coolant flow rate measurements, 
d.) rotational speed m.easurements , and e.) electrical power measurements. 
A schematic drawing of the Instrumentation setup Is shown In FIGURE 2. 
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a,) Temperature Measurements 

Temperature readings are obtained from thermocouples 
mounted in the evaporator wall to determine a temperature profile 
and associated heat flux through the wall. Pour grounded copper- 
constantan thermocouples with a l/l6-inch diameter stainless-steel 
sheath are utilized. They are placed within l/l6-inch diameter 
wells drilled 5j-inches deep into the stainless-steel evaporator. 

To insure that no minute air gaps are present between the well walls 
and the sheaths, a coating of Silver Goop is placed on the thermo- 
couple prior to insertion. Two of the thermocouples are mounted 
opposite each other at the same radial position to check on the 
uniformity of the temperature field. This pair is mounted at the 
middle of the ^-'*nch wall thickness. The other two thermocouples 
arc also mounted opposite one another, one l/l6-inch radially from 
the outside evaporator suface, one l/l6-inch radially from the inside 
surface , 

Wie leads from the thermocouples are brought through both 
evaporator and forward condenser flanges, down the side of the condenser, 
and through the rear condenser flange to phenolic Junction boards. 

From the boards the leads go inside the shaft and down to the slip-ring 
unit which is attached to the main assembly by a coupling mounted at 
the end of the shaft. 

A thermocouple is suspended in the liquid annulus of the evaporator 
to measure the bulk liquid temperature during operation. A sixth 
thermocouple is suspended in the vapor space at the evaporator exit to 
record the saturation temperature of the vapor. The leads to the 
thermocouples in the liquid and vapor pass through slots in the evaporator 
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wall (see photograph in FIGURE 4) after passing through the condenser 
forward flange. 

The temperature of the coolant is measured before and after spraying 
by quartz themoraeters mounted in mixing tubes located above and below 
the condenser spray cooling cylinder. 

b. ) Pressure Measurements 

The saturation pressure of the vapor is measured by a miniature 
semi-conductor pressure transducer mounted in the vapor space. The 
transducer is temperature compensated and is threaded on a ^-inch 
diameter arm that extends from the shaft through the end plug and 
into the vapor space. The probe's leads pass through the arm, up 
through the shaft to the junction boards, back into the shaft, and 
down to the slip-rings. 

The pressure of the coolant prior to spraying is indicated by 
a pressure tap on the coolant feed line. 

c. ) Coolant Flow Rate Measurements 

The coolant flow rate is obtained by passing the coolant 
through a flow rotameter prior to spraying. 

d. ) Rotational Speed Measurements 

The rotational speed of the heat pipe is found by two different 
methods. First a Hewlett-Packard optical tachometer is employed. Six 
equally spaced spots mounted circumferentially on the shaft cause the 
tachometer to emit an electrical impulse each time one passes through 
its optical field. These in^julses are fed to a Donner frequency 
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counter and counted over a ten second interval. The read-out Is in 
revolutions per cdnute. The second method of determining the rpm 
is by the use of an electronic strobe light, 
e.) Electrical Power Measurements 

The amount of electrical power that the heater is using 
is found by placing a voltmeter in series and an ammeter in parallel 
with the power lines returning from the heater. 

The power level is controlled by a field rheostat placed 
in parallel with the DC motor-generator. A 150-ampere circuit breaker 
is Included in the circuit, 

A slip-ring unit must be used to obtain readouts from the 
instrumentation mounted on the rotating heat pipe. A photograph 
of the slip-ring unit is shown in FIGURE 6- It utilizes liquid mercury 
as a viscous contact medium between the revolving and stationary leads. 

This cuts the frictional heating noise to less than 10 microvolts at 
speeds up to 4000 rpm. The slip ring unit has 22 terminals and is 
therefore capable of accomodating 9 thermocouples and 1 pressure 
transducer. The output from the slip-rings is fed to a Hewlett-Packard 
2010 C Data Acquisition System containing an integrating digital 
voltmeter, guarded data amplifier, and guarded cross-bar scanner. This 
system provides printed digital output with an accuracy of + 0.5 microvolts. 

C. PROPOSED OPERATING PROCEDURES 

The operating procedures are divided into three sub-procedures: 

1.) filling procedure, 2.) testing procedure, and 3.) treatment of 
data. 
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1 . ) Filling Procedure 



Before operating and testing the rotating heat pipe, the 
working fluid must be Introduced Into the system. This Is done In 
such a manner as to keep the amount of non-condensable gases In the 
system to a minimum. These gases reduce the amount of condenser 
surface available to the vapor. Before filling commences, a maximum 
operating saturation pressure Is decided upon. A limit on the operating 
pressure Is Imposed by the amount of pressure the heat pipe Is designed 
to withstand. The maximum operating pressure determines the minimum 
amount of working fluid that Is added to the pipe (the maximum Is 
determined by the amount of liquid the evaporator annulus Is capable 
of holding) . 

The pipe is placed In a vertical position with the evaporator 
section down and the shaft not In place. The working fluid is added 
and a heat flux is generated by passing a current through the heater 
element. As the liquid evaporates, the vapor drives the non-condensable 
gases to the end of the condenser and out through the hole In the end 
plug. Once It Is felt that the gases are eliminated from the system, 
the shaft Is screwed into the threaded shaftway In the end plug until 
it seals against a metallic 0-ring located on its end. Heat Is still 
added and liquid evaporated to determine If there is enough working 
fluid to attain the pre-determined maximum operating saturation pressure. 
Once this is ascertained, the current to the heater Is shut off. The 
filling phase of operation is completed and testing can commence. 

2 , ) Testing Procedure 

Once filling is completed, the test stand is placed in the 
chosen orientation for the run. Rotation Is started and the rpm is 
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brought to a pre-deterrained level for the run. This is checked by 
readings obtained from the optical tachometer and strobe light. 

The spray cooling system is turned on. Finally, a current is passed 
through the heater element. The current is increased in steps with 
a waiting period between each increase to allow the pipe to reach steady- 
state. The readings of the thermocouples, ammeter, voltmeter, quartz 
thermometers, pressure transducer, coolant pressure tap, and rotameter 
are recorded with each current step. This is continued until the 
maximum saturation pressure level is reached. The same procedure is 
then repeated for a series of increasing rpm values. 

Shut-down is accomplished by first shutting off the current to 
the heater, then decreasing the rpm by use of the remote control unit 
(simply turning the motor off will engage the magnetic disc brake). 
Finally, the coolant system is secured. 

3 . ) Treatment of Data 

The amount of heat the heat pipe is transporting during 
steady-state operation is the primary quantity to be determined from 
the experimental data. Because of its importance in fixing the 
operating characteristics of the rotating heat pipe, three methods are 
used to obtain the heat flux. 

The first method is to determine a radial temperature profile 
in the evaporator wall from the readings of the thermocouples imbedded 
in the wall. The heat flux is then obtained from Fourier's heat 
transfer equation in cylindrical coordinates. 

The heat flux is also obtained by recording the voltage and 
current levels to the heater. Once the boiler heat losses are 
approximated, either by the analytical approach described in Appendix 
A (using different convection coefficients calculated for different 
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speed levels) or the experimental procedure described in ref. [6] , 

the actual heat flux is obtained by subtracting this loss from the 
gross power input. 

The third method for determining the heat flux is by performing 
an energy balance on the condenser cooling section. The coolant 
flow rate is obtained from the flow rotameters and from a calibrated 
curve for the standard nozzles spraying at a fixed pressure, known 
from the pressure tap reading. The coolant's mean temperature is 
measured before and after spraying by quartz thermometers mounted 
in mixing tubes located above and below the spray cooling cylinder. 

With this information and knowing the coolant's density and specific 
heat, a heat flux is obtained from a simple energy balance. 

The vapor pressure is another important quantity that is 
monitored during heat pipe operation. It acts as the independent 
variable in the experimental analysis of the rotating heat pipe. The 
vapor pressure is determined both from the reading of the pressure 
transducer in the vapor and the thermocouple mounted in the vapor 
space. The thermocouple gives an indication of the saturation 
temperature of the vapor, which can be converted to an equivalent 
saturation pressure. 

The rotational speed is the adjustable parameter in the 
experimental study of the rotating heat pipe. As indicated previously, 
its value is obtained by both an optical tachometer and electronic 
strobe light. 
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III. ANALYTICAL PROGRAM 



A. PRELIMINARY ANALYSIS 

In his thesis on the operation of a rotating wickless heat pipe, 
Eallback J7J discussed the limitations on the rotating pipe due to 
various fluid dynamic phenomena. Using experimental correlations 
taken from systems approximating the dynamics of the new heat pipe 
he was able to estimate the Timits of operation imposed by the mechanisms 
of entrainment, sonic velocity, and critical nucleate boiling. In 
order to obtain a condenser limit for the conically-shaped condenser 
section Eal?back performed a Nussel t-tji)e analysis for the condensa- 
tion on a truncated cone. Existing solutions for rotating condensation 
[8, 9J were not applicable to this geometry since they could not satisfy 
the condition of zero condensate velocity at the truncated end of the 
condenser. Eallback 's basic assumptions for such an anlysis were 
justified by classical theory or experimental correlations. With these 
assumptions and the coordinate system shown in FIGURE 8, he performed 
force balances in the X-and Y-coordlnate directions on an infinitesimal 
fluid element. After integration and simplification of these basic 
momentum relationships, an expression for the condensate velocity was 
obtained . 

2 T 

u„ = ^ y - [6Y - ^ ] [(R+Xsin(j-6 cosCJ) (slnd-cos0 dX )] (1) 

where u = fluid velocity, (ft/hr), 

= fluid density (Ibm/ft^) 

= fluid dynamic viscosity, (?bm/ft-hr) 

CO = rotational speed (l/hr) 

Y = coordinate measuring distance normal to surface, (ft) 
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5 = condensate film thickness, (ft) 

= minimum wall radius in the condenser section, (ft) 



0 = half-cone angle, (rad) 

X = coordinate measuring distance along the surface, (ft) 




slope of the condensate film 




FIGURE 8 Coordinate System for Condensing Limit 
From continuity^ an expression was obtained for the condensate mass 
flow rate. 

Ballback then assumed that the film is very this (ie,5«R,) and that 

(3) 

This assumption allowed the mass flow rate expression, eq'n(2), to 
be simplified to 

I [(fi+Xs.njz()^,n^] (4) 



With a steady-state assumption, the rate of condensation equation 
became 



dX '■<3X 




St, . 
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(5) 



where 



h 



fg 



= latent heat of vaporization of the working fluid 
(Btu/lbm) . 



/ 



■^t was obtained by considering the heat transfer by conduction 
dX 

to an element of surface area. In determining this expression, 

Ballback considered only the thermal resistance due to the condensate 
film and neglected the thermal resistances in the condenser wall and 
outside cooling mechanism. This simplified expression was substituted 
in equation (5) which after integration, yielded a closed-form solution 
for the condensate film thickness as a function of X. Ballback defined 
a local heat transfer coefficient as the thermal conductivity of the 
condensate divided by its thickness. He integrated this expression 
over the surface of the condenser to obtain a mean heat transfer 
coefficient. With this mean coefficient a total heat rate for the 
condenser was derived. Ballback then coupled this condenser heat rate 
with a conduction equation across the condenser wall to generate a 
condensing limit curve. 

Ballback 's results were obtained by neglecting the themal 

» 

resistances in the condenser wall and the external condenser cooling 

mechanismj and by making the additional simplifying' assumption that 

^<< tan^ . (3) 

dX 

It is apparent that the above assumption, eq'n(5)} is true only for 
large half-cone angles. Also if is zero the mean heat transfer 
coefficient obtained in Ballback' s analysis is identical to coefficient 
obtained in an existing boundary layer analysis of condensation on 
a rotating cone performed by Sparrow and Hartnett [9 ] . They point out 
that their solution is limited to large half-cone angles. 
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Ballback's theory can be improved if it is modified to take 
into account the thermal resistances in the condenser wall and 
condenser cooling mechanism. This is accomplished by including these 
teirnis in the incremental heat conduction equation, which, for a thin 
condensate film and condenser wall, becomes 



dQ^ 2 7T (R+Xsind) (T -T, J 
-T^t= o ' s inf 

i»/k +t/k + 1/h 
f w 



( 6 ) 



where k = thermal conductivity of the liquid condensate 
(Btu/hr-ft-°P) 

k = thermal conductivity of the stainless-steel 
condenser wall (Btu/hr-ft-°F) 
t = thickness of condenser wall, (ft) 
h = condenser cooling mechanism convection 
coefficient, (Btu/hr-ft -°F) 

When this expression is substituted in the rate of condensation equation 
(5) along with Ballback's simplified expression for the condensate mass 
flow rate, a first-order, non-linear differential equation for the 
condensate film thickness if generated. A closed-form analytical 
solution for 6 cannot be obtained from this formulation. A Runge-Kutta- 
Gill numerical Integration scheme was therefore used to solve eq'ns(4), 
(5) and (6) using an IBM - 3^0 Model 67 digital computer and Ballback's 
heat pipe dimensions. A zero initial condensate film thickness was 
assumed to start the integration scheme since this is the only value 
that gives a zero condensate velocity at the truncated end of the 
condenser. Since this theory is still limited to large half-cone 
angles by the neglect of the condensate film slope in the mass flow 
rate expression^ the results obtained from the above numerical 
integration are approximate and are assumed to represent an upper limit 
on the condensing limit for rotating heat pipe operation. 
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The condenser limit heat transfer rates are plotted versus 
rotational speed for a saturation temperature of 212^F and half- 
cone angles of one— and two— degrees in FIGURES 9 and 10. The 
outside convection coefficient is a parameter showing the 
dependence of the heat transfer on the condenser surface cooling 
mechanism. With an infinite convection coefficient and a 3 half- 
cone angle, a comparison between Ballback's predicted condensing 
limit and the condensing limit calculated from the modified theory 
is shown in FIGURE 11. The modified theory is seen to predict a 
slightly lower condenser limit heat transfer rate; apparently, 
neglecting the wall resistance during integration leads to higher 
rates of heat transfer. 

B. REVISED ANALYSIS 

The large half-cone angle restriction on the above modified 
analysis can be removed by including the film thickness slope term 
in the expression for the condensate mass flow rate. 



When this expression is differentiated and substituted in the rate of 
condensation equation for steady-state operation eq*n(5), a second- 
order differential equation for the condensate film thickness is 
obtained . 




(a) 
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[r+Xs.v,<J^Scc 4] [s»>f C»4^] jl'^'t']] 



2TT( X^l j>) -Tiyfi) 

Vk^ + t/Kw+ '/h 



( 7 ) 



This expression is valid for all half-cone angles. It reduces upon 
simplification to existing theoretical expressions for a rotating 
disc [9] at 0=90^ and for a rotating drum [10] at 0=0^. 

A numerical integration method must be utilized to solve 
equations (6) and (7) simultaneously to obtain a solution for the 
condensate film thickness as a function of X and the condenser limit 
heat transfer rate Q^. To start the integration, initial values for 
the heat rate, film thickness, and film thickness slope must be 
determined. is zero at X=0 since there is no area for heat transfer 

to occur (the truncated base is assumed to be adiabatic) . Leppert 
and Nimmo [10] state that for condensation on finite surfaces normal 
to inertial forces, the starting condensate film thickness is a 
function of the minimum film thickness, which occurs at (or very near) 
the surface’s edge. This is shown in FIGURE 12, 




FIGURE 12 Behaviour of Condensate Film Under Influence of Inertial 
Forces 
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The minimum thickness depends on the particular overfall condition 
occurring at the surface's edge. Clifton and Chapman[12] use a 
free-overfall edge condition to analyze natural convection on a 
finite horizontal plate. If density variations are replaced by 
rotational acceleration as a driving force, an expression can be 
obtained for the minimum condensate film thickness using a free- 
overfall edge condition. This analysis is outlined in Appendix C 
with the result that 



equation (1) in order for the condensate velocity to be zero at 

the truncated end of the condenser section, the slope of the condensate 

film must initially equal the tangent of the half-cone angle. 



The method of solution of the differential equations for the 
condensate film thickness and condensing limit heat rate resembles 
that for the solution of a system of non-linear differential equations 
with an asymptotic boundary condition. With the initial values of 
the heat rate and condensate slope known, an initial value of the 
condensate thickness is guessed to start the numerical integration 
scheme. With these starting values equation (7) is integrated, 
determining the condensate film thickness as a function of X. The 
film thickness at X=L^ is compared to that obtained from the free- 
overfall condition eq'n(8). The overfall thickness is a function 

iBecause of the coordinate system Csee FIGURE 8) , this is an approx- 
imation for all half-cone angles but zero. The exact boundary condition 
for these angles is not known because the condensate velocity is not 
zero at X=0 




.where S =condensate film thickness at the overfall, X=L . From 
c c 




(9) 
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of the condensate mass floi-r rate at X-L , which is calculated from 

c 

equation (2) . If the two tliicknesses compare, a solution is 
obtained. If they do not compare, a new guess for the initial film 
thickness is used as a starting value in the integration scheme. 

This procedure is repeated until a satisfactory result is arrived at. 

A Runge-Kutta-Gill numerical integration scheme was used to 
solve equations (6) and (7) simultaneously on an IBM 360 Mod 67 
digital computer. The initial condensate film thickness was iterated 
until the final thickness obtained from the integration and that 
obtained from the overfall condition agreed to within .0004 inches. 

Very favorable results were obtained from this procedure for 
a half-cone angle of zero degrees. However, as the value of the 
half-cone angle was increased, a problem was encountered which 
centered about the instability of the highly non-linear second-order 
differential equation for the condensate film thickness. As is the 
case with more well-known non-linear differential equations with 
asymptotic boundary conditions [14] the instability of the equation 
was highly sensitive to the starting value for the film thickness. 
Consequently, no numerical results were obtained for half-cone angles 
greater than zero. The condensing limit heat transfer rate is 
plotted versus rotational speed for a zero half-cone angle for different 
values of the outside condenser surface convection coefficient, in 
FIGURES 13, 14, and 15. The saturation pressure and therefore temperature 
is a parameter on these graphs. When the thermal resistances in the 
wall and outside cooling mechanism are neglected, the results compare 
with those of Leppert and Nimmo [10] for a rotating cylinder. 
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IV SUMMARY 



A . EXPERIMENTAL GEOMETRY 

1. The design and machining of the rotating heat pipe has been 
accomplished as described earlier. It is designed such that visual 
observations can be made of internal pipe mechanisms such as critical 
nucleate boiling and condensate entrainment. The instrumentation system 
is devised so that the parameters of the experimental study such as 
heat flux^ saturation pressure, and rotational speed can be monitored 

in more than one manner so that checks on a reading can be obtained. 

A high degree of precision was used in the machining of the pipe 
and the mounting of the instrumentation to keep the degree of un- 
certainity in the determination of the operating characteristics to a 
minimum , 

2. The heat pipe and drive assembly are mounted on a special test 
stand so that tests can be taken on the device at different orientations 
in the gravity field. 

3. Safety is insured during operation by enclosing the rotating 
assembly within stainless-steel shielding, Plexi-glass viewports and 
doors in this shielding still allow visual observation of the operation 
to be accomplished, 

B. THEORETICAL STUDY 

1. The preliminary theoretical investigation performed by Ballback 

is limited to large half-cone angles and by the fact that the mean 
heat transfer coefficient was obtained by considering the effect of 
the thermal resistance due to only the condensate film and neglecting 
the effects of the condenser wall and outside cooling mechanism. 
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2, After amending Ballback's theoretical results to include the 
thermal resistance effects of the condenser v;all and outside cooling 
mechanism, the condenser limit for a particular operating condition is 
seen to decrease (i.e., less heat can be removed than was previously 
expected) , 

3. Based on Ballback's basic Nusselt-type analysis of condensation 
on a truncated cone under the influence of rotational forces, an 
improved theory was obtained that is good for all half-cone angles and 
includes the effects of thermal resistance in the condenser wall and 
outside cooling mechanism. This formulation was coupled with a free 
overfall condition and produced results for a 0° cone angle that are 
in close agreement with those from existing theory for a rotating 
cylinder. Results were not obtained for non zero half-cone angles 
because of the instability of the governing equation, the dependency of 
this condition being on the accuracy of the Initial guess for the film 
thickness at the truncated end of the condenser section, 

4, The condenser limit of the rotating heat pipe is govemed by 
the maximum pressure the pipe is designed to withstand and the amount 
of working fluid it can contain. If this peak operating level is 
designed high enough the performance of the pipe may become limited by 
entrainment or critical nucleate boiling before the condenser limit is 
approached . 

5. The performance of the rotating heat pipe is governed by the method 
by which the condenser section is cboled. 
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V. RECOMf^lENDATIONS 



A. EXPERIMENTAL ASPECTS 

1. The heat pipe assembly should be completed and the heat transfer 
capabilities tested for various operating conditions as described in 
this study, 

2. Different condenser geometries should be designed, machined, and 
tested. Some of the variations could Include the varying of the half- 
cone angle, the addition of some non-wetting material to the condensation 
surface to promote dropwise rather than film condensation, the machining 
of longitudinal grooves on the condensing surface to decrease the effect 
of the interfacial shear on returning condensate, and the placing of 
fins on the outside condenser surface to increase the effective heat 
transfer area, 

3 . Different condenser cooling mechanisms, from air to liquid gases 
should be utilized and studied, 

B. THEORETICAL ASPECTS 

1, An iterative scheme such as that described ref, [141 should be 
formulated and included in the conputer program to aid in obtaining 
theoretical results for half-cone angles other than zero degrees. 

2, The effect of the vapor should be included in the theoretical 
analysis by solving the continuity, momentum and energy equations in 
the vapor space and coupling this solution with that of the condensate 
by the common shear stress and temperature at their Interface. 



49 



APPENDIX A 



DETERMINATION OF MAXIMM RADIAL 
HEAT FLUX AND ASSOCIATED LOSSES 



1, The heater element is made of Chromel A which has a melting 

point of approximately 2400°F. It is assumed that its maximum 

0 

surface temperature is 2000 F. Using a cylindrical coordinate system 
I = 2 WA T/R^ 

\‘^^ss'^^inc‘''^0 

where 



R ^thermal resistance of stainless steel 
boiler wall (°F-hr-f t/Btu) 

R, •--•thermal resistance of inconel heater 
^^^sheath (°F-hr-f t/Btu) 

R ^thermal resistance of magnesia heater 
^insulation, (°F-hr-f t/Btu) 

The resistance of the fluid annulus is neglected due to boiling. From 
the known geometry of the evaporator section and the known thermal con- 
ductivities of the various elements 

1,8125 °F-hr-ft 

R = 1,5625 =.01237 Btu 
ss ^2 

, 1,8325 °F-hr-ft 

R, =1,8125=. 00127 Btu 

^ ^1.8 32 5 °F-hr-ft 

Rj^0= 1-8325 =.01275 

1,2 

R^=,02638°F-hr-ft 
^ Btu 

Assuming a maximum operating saturation pressure of 3 atmospheres, the 
temperature difference across the resistance is 
T=2000-270=1730°F 
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Then 



Btu 

hr-ft 



Q _ 2 TT(1730) 
L “ .02638 



420,000 



The heating surface Is approximately 4 inches long and the boiler inside 



surface area is 

A^=2 (1.5625) (4) (1/144)=,275 ft ^ 



The heat flux is then 

q= j ^ 470 , 000 



Btu 

hr-ft' 



2. The radial heat losses are calculated in much the same manner. 



Neglecting the effect of the thin coating of Sauereisen cement, the 



thermal resistance is 



R,“R, +R,, ,.,+R^ + 

t inc MgO ins 



1 

hA 



where 

R^^^=thermal resistance of insulation, (°P-hr-ft/Btu) 

h=outside convection coefficient between^ the 
insulation and atmosphere (Btu/hr- P-ft ) 

A=outside surface area of insulation. 

It is assumed that the heat transfer on the outside of the 
rotating evaporator assembly can be approximated by the flow over 
a stationary cyliner. Then ref. (11 ] can be used to obtain an 
estimate for the outside heat transfer coefficient. 



NUp=C(Rejj)" 

Rgp= 200,000 

c= .0239 
m= .805 

NUp=362 

h =79 Btu/hr- °P-ft^ 

.0277 °P-ft-hr/Btu 

2.00 

R, = ^'^1.98=. 001 16 °P-hr-ft/Btu 
inc — 75-^ ' 



for w =3600 rpm 

Rejj =(3600) (27 t) ( 60) (d) 

V . 
air 

where d=diameter of insulation, 

(ft) 

V . =kingmatic viscosity of air 
(ftVhr) 




1.980 

1.952 = .01197 
TT2 



°P-hr-ft/Btu 
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1 



w 



1 2.75 

„= 2.00 =6.37 °F-hr-ft/Btu 

ins 

These v/ere obtained by assuming a 3/4-inch layer of insulation with 
a thermal conductivity of .05 Btu/hr-°P-f t^ 

R^=6.41 °F-hr-ft/Btu 
I ^ ^ ' ^ ' 6^4 ' f '" Btu/hr-ft 
’loss"! ' ”50 Btu/hr-ft" 

h can be re-calculated for different rotational speeds, giving a 
relationship between heat loss and rpm. 

3. Using the known evaporator geometry and assuming that I9 helical 
turns of heater can be wrapped on the 4-lnch heating surface, the total 
length of heater is 

2 IT (1.90625) (19)/12= 19 ft 

The resistivity of the heater element is known to be .0785 ohms/ft. 

This fixes the total resistance in the circuit 
R=( .0785) (19)=1. 4887 ohms 

If the maximum voltage that can be supplied is 220 volts, the current 
is 

1=220/1.4887= 150 amperes 

The gross power input is therefore 
P=I^R=32.5 kilowatts 

2 

which corresponds to a heat flux of approximately 350; 000 Btu/hr-ft , 
assuming losses of lO/. 

4. It is therefore seen that the total heat flux capable of being 
conducted across the evaporator wall is limited by the amount of heat 
that can generated by the heater. 
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APPENDIX B 



PRELIMINARY UNCERTAINTY ANALYSIS 



Assume that the heat flux is to be determined at a surface at 

a radial distance R by the readings of two thermocouples placed at 

radial positions r^ and • The governing equation is 

o 1 



9 -JiiZallil 

“I ° Ar ■ R In Cn.'/r,') 



-where 



T 

0 

T. 

1 






k = 



“ temperature reading obtained at r^ 

= temperature reading obtained at r^ 

= surface area at the radial position R 
thermal conductivity of material being used 



Then 








"From the governing equation 




hq , ^ 





-Ka-Tf ) 

WK/r{)] 



X 






« ‘ R* lvi(n;A'i} 

rj Rpii W/n;^ 



-Then it is seen 

^ ^Tf I 

R ■ K fqpir) ■ a-T,) * R ■ r; lti(Vr,') r'l 

If it is as sun -d that all the uncertainties are additive and that the 
values of k and R are known with high accuracy, it is seen that the 
uncertainty in the determination of the heat flux depends on the 
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uncertainty in the readings of the thermocouples and in their positioning. 
It is also seen that the farther apart the thermocouples are placed, 
the larger the value of the temperature difference and the logarithmic 
terms become and the smaller the effect of their associated uncertainties 
become. 
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APPENDIX C 



DERIVATION OP FREE 
OVERFALL END CONDITION 

From ref. [l3]it is seen that during flow down a mild, inclined slope, 
critical flow will occur at a sudden change in slope as shown below. 




From energy considerations of the system 



dx 



(1-Pr^) 



S - 
0 



S 



f 



where 



Fr = Froude Number = 
S 



V 
/ET 
dz 



bed slope 



= friction slope = 



c^R 



R = wetted mean radius (=A/P) 
s 

c = Chezy constant (similar to friction factor, f) 

In the region of mild slope, since the slope is so mild. On the 

other hand, in the region of steep slope, S^<S^. Therefore there 
must exist a point, or region, where . Here 

4? (1-Fr^) = 0 
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Since the slope of the free surface is not zero, it follows that at 
this point 

Pr = 1.0 

which is the critical condition. Thus, 
v/yiF^ = i.o 
or 

= vVg 

where is the critical depth. 

For flow over a sharp corner^ the flow can separate from the 
steep wall as shown below. 




/ 

/ 



This is called a free overfall. In this case^ there is experimental 
data [15 ] to show that 

The value of occurs very near the edge of the plate_, so that in the 
case of the condensate film flowing on the condenser wall 

ati = L^ 
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where 



=5 condensate film thickness at the overfall 
« the condenser length 
7c " V^/g 

Here, V is the average velocity of the condensate given by 




where 

= velocity of the condensate 
= condensate mass flow rate 
Pf = density of the condensate 

Also, the acceleration of gravity, g, is replaced by the rotational 



acceleration term 

g = <0 ^ R 

where 



CO = rotational speed 

= small radius of the condenser cone 
0 = half-cone angle 



Therefore, 



or j at X=L 



yt\f 

^ sin<^> ) (2 itRS 



0.7 rrif 









Finally, as the free overfall condition in the condenser. 




0-7 M 

f2TT)^ (R^+L^sin</))'^ 
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